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ABSTRACT 

It may be possible to construct a laser interferometer gravitational wave antenna in space with hmis ^ 
10~23 at / - O.lHz in - 2020. This deci hertz antenna may be called DECIGO/BBO, which stand 
for DECi hertz Interferometer Gravitational wave Observatory and Big Bang Observer, respectively. 
The analysis of 1-10 years observational data of the coalescing binary neutron stars or black holes at 
the distance of SOOMpc will give us the spatial position within ^ an arc minute and the time of 
the coalescence within ^ 0.1 sec beforehand. With the knowledge of the accurate position and the 
time of final merging event, the follow up simultaneous observation using high frequency (/ lOOHz) 
gravitational wave antennae as well as electro-magnetic wave antennae from the radio frequency to 
the ultra high energy gamma ray will reveal the physics in the enigmatic event of the coalescence and 
the formation of the black hole. 

Subject headings: binaries:general - gravitation - gravitational waves 



1. INTRODUCTION 

There are at least four frequency bands where active 
or planned gravitational wave antennae exist. The reso- 
nant bars as well as ground laser interferometers such as 
TAMA300 , LIGO I , VIRGO , and GEO600 are covering 
the frequency band of lOHz-kHz while the timing analy- 
sis of pulsars is covering ~ lO^^Hz band. Laser Interfer- 
ometer Space Antenna ^ (LISA) will cover lO"** — 10~^Hz 
from ~ 2010. Very recently the space antenna in the deci 
Hertz band (10~^ — lOHz) comes into the group. Accord- 
ing to Seto, Kawamura and Nakamura (2001) it might be 
possible to construct a laser interferometer gravitational 
wave antenna in space with hrms ^ lO^^*" at / ^ O.lHz in 
this century. Using this antenna they show that 1) ~ 10^ 
chirp signals of coalescing binary neutron stars per year 
may be detected with S/N ~ 10**. 2) The time variation 
of the Hubble parameter of our universe may be deter- 
mined for ten years observation of binary neutron stars 
at z 1 so that we can directly measure the accelera- 
tion of the universe. 3) the stochastic gravitational wave 
which is predicted by the inflationary universe paradigm, 
may be detected. 

Seto et al. (2001) call the deci Hertz antenna DE- 
CIGO(DECi hertz Interferometer Gravitational wave 
Observatory) while recently in the NASA SEU 2003 
Roadmap "Beyond Einstein" the deci Hertz antenna is 
called as BBO (Big Bang Observer) ^ to stress the detec- 
tion of the stochastic gravitational waves from inflation. 
Since DECIGO and BBO are similar antennae, in this 
Letter we call the deci hertz antenna as DECIGO/BBO. 
In DECIGO, hrms ^ 10"^^ around 0.1 Hz is the ulti- 
mate goal in the far future, assuming the quantum limit 
sensitivity for a 100 kg mass. 

Let us consider more realistic parameters which might 
be achieved in ~ 2020 such as 300W laser power, 3.5m 
mirrors, 5 x lO'^km arm length and 0.01 LISA accelera- 
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tion noise, which will give hrms ^ 10~^^ around 0.1 Hz 
(e.g. Larson, Hiscock & Hellings 2000). In this sensitiv- 
ity, S/N of the chirp signal from the coalescing binary 
neutron stars at z ~ 1 will be ~ 1 while that of IOMq 
black hole binary will be ~ 10 so that we may directly 
measure the acceleration of the universe similar to the 
ultimate DECIGO. After subtracting the signals from 
neutron star binaries and black hole binaries ^ we may 
detect the stochastic gravitational waves from inflation 
if ^Iqw ^ 10"^^ which is the upper bound from CMB 
quadrupole anisotropics measured by COBE (Gorski et 
al. 1996). 

One of the main targets of the ground laser interfer- 
ometer gravitational wave antennae is the chirp signal 
from the coalescing binary neutron stars or black holes 
at the distance of 200 — 400Mpc with the expected event 
rate of ~1 per year (e.g. Phinney 1991). Since the event 
rate is ^^1 per year, the same binary is in the band of 
DECIGO/BBO, 1 10 years before the ground laser in- 
terferometers and bar detectors detect the gravitational 
waves from the final merging event. This means that in 
practice we may observe the evolution of the coalescing 
binary from 10^^ Hz band to ~kHz band. The duration 
in the band of the ground laser interferometer is three 
minutes or so while the duration in DECIGO/BBO band 
is 1-10 years so that we may extract the information of 
the binary from the DECIGO/BBO band observation 
first. The point here is that we may know the accurate 
position and the distance to the binary beforehand so 
that we can prepare various detectors for the observa- 
tion of the final merging phase, that is, we can point 
electro-magnetic antennae from radio to ultra high en- 
ergy gamma rays and cosmic rays to the direction of the 
coalescing binary at the expected final merging time. We 
may also tune the gravitational wave antennae to vari- 
ous sensitivity bands to detect the characteristics of the 
merging event such as the frequency at the inner most 

^ There may be other sources that we do not know about at this 
point. 
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stable circular orbit and the quasi-normal mode of the 
final black hole. 

In this Letter using the practical DECIGO/BBO in 
~ 2020 we show how accurately we can determine the 
spatial position and the distance to the coalescing binary 
neutron stars or black holes. Since the angular resolution 
of LISA was calculated in (Cutler 1998; Hughes 2002; 
Seto 2002; Vecchio 2003), we apply these methods to the 
case of DECIGO/BBO. We adopt the Hubble parameter 
h = iJo/lOOkm/s/ Mpc=0.7 and the units of c = G = 1. 

2. GRAVITATIONAL WAVEFORM AND PARAMETER 
ESTIMATION 

2.1. Sensitivity of DECIGO/BBO 

DECIGO/BBO would consist of three spacecrafts sep- 
arated by 5 X 10'* km which is 1/100 times smaller 
than the size of LISA (see Bender et al. 2000). They 
would be spaced in an equilateral triangle and be orbit- 
ing around the Sun. The change of the detector's orien- 
tation and the position enables us to obtain the source 
position in the following two ways; (i) DECIGO/BBO's 
orientation rotates with a period of one-year, which im- 
poses the modulation on the measured signal, (ii) DE- 
CIGO/BBO's orbital motion around the Sun imposes the 
periodic Doppler shift on the signal frequency (see also 
Cutler 1998 in the case of LISA). For the monochromatic 
signal of the frequency /, (i) the rotational modulation 
changes the frequency / to / ± 2/r, where T = 1 yr 
(2/T ~ 10~^ Hz), (ii) the Doppler modulation changes 
/ to /(I ± ti), where v ^ 10^^. Then for the frequency 
higher than ~ lO'^ Hz like in DECIGO/BBO, the effect 
of the Doppler modulation is more important for the de- 
termination of the source position than the rotational 
modulation so that we only consider the former effect in 
the measured signal. The observed waveform h in the 
frequency domain h]j(f) is given by. 



in the frequency domain is given by. 



hD{f) = h{f) e^*-^f^ 



(1) 



where h{f) is the waveform at the solar system barycen- 
ter. (poif) is the Doppler phase, which is the differ- 
ence of the phase between the detector and the Sun; 
0^ = 27r/i?sin6'scos(27rt/r - (j^s), where R = I AU, 
T = 1 yr and {Os,4's) is the direction to the source. 
These angular coordinates are defined in a fixed barycen- 
ter frame of the solar system (see Cutler 1998). 
The strain sensitivity of DECIGO/BBO is about 1000 
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times better than that of LISA (i.e. 10 Hz 
/ = 0.1 — 1 Hz), and the acceleration noise is 100 times 
lower than that of LISA (i.e. 3 x lO-i^m/s^/Hz^/^). 
Adopting these parameters, we show the sensitivity of 
DECIGO/BBO as compared with LISA and LIGO II in 
Fig.l (Larson et al. 2000). Note that we are using a 
smooth curve and neglecting the wavy behavior of the 
transfer function for simplicity. 

2.2. Gravitational Waveform 

We consider the equal mass neutron star and black hole 
binaries as sources of DECIGO/BBO. We only consider 
a circular orbit since the expected eccentricity is ~ 10"'^ 
for the neutron star binary. We use the restricted 1.5 
post- Newtonian approximation as the in-spiral waveform 
(Cutler & Flanagan 1994) for simplicity. The waveform 
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where A is the amplitude and ^'(/) is the phase. They 
depend on six parameters; the red-shifted chirp mass 
M, = (MiM2)3/5(Mi + M2)-^/^{l + zs), the reduced 
mass fiz — MiM2{l + zs)/{Mi + M2), the spin-orbit cou- 
pling constant the coalescence time tc and the phase 
(f)c and the luminosity distance to the source Ds- The 

amplitude is given by ^ = MI^^ /{tt^/^Ds), 

where K is the constant determined by the inclination 
of the source, the relative orientation of the source and 
the detector. Since the average value of K is about unity 
(Finn & Chernoff 1993), we assume K = 1 for the fol- 
lowing calculation. The phase ^(/) in Eq.® is a rather 
complicated function of A^^, fi^, and tc which is 

given in Eq.(3.24) of Cutler & Flanagan (1994). 

2.3. Parameter Estimation 
The signal observed by DECIGO/BBO, hoif), can be 
obtained inserting Eq.(|2l to Eq.QJ. The signal hnif) is 
characterized by eight parameters {Aiz, ^^z, /?, i'c^tc, Ds, 
6st 4>s)- In the matched filter analysis with the template, 
these parameters can be determined. We compute the 
errors in the estimation of these parameters using the 
Fisher matrix formalism (Finn 1992; Cutler & Flanagan 
1994). The variance-covariance matrix of the parameter 
estimation error is given by the inverse of the Fisher 
information matrix Tij as (A7,;A7j) — (r~^) 
Fisher matrix becomes 
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df dhUf) Ohoif) 



(3) 



Sn{f) dji djj 

where Sn{f) is the noise spectrum. We regard Sn{f) 
as the instrumental noise in Fig.l, neglecting the binary 
confusion noise since there are no or little confusion noise 
for / > 0.1 Hz (Seto et al. 2001). The signal to noise 
ratio {S/N) is given by 



(S/N) 



df 



Snif) 



hoif) 



(4) 



We integrate the gravitational waveform in Eq.® and 
Eq.|0J) from 1 or 10 yr before the final merging to the cut- 
off frequency fcut when the binary separation becomes 
r ^ 6(Afi + M2). 

The initial frequency is given by /mit = 
0.23(MjMQ)-5/8(To6,/lyr)-3/8Hz where T^bs = 1 
or 10 yr, and the cut-off frequency is fcut — 
4.4 X 10^[{Miz + M2z)/M0]-iHz. The resuU does 
not depend on the value of /c^t much, since Sn(f) is 
large at /c^f 

3. RESULTS 

We consider the neutron star binaries (I.4M0-I-I.4M0) 
and the stellar mass BH binaries (IOM0 -I- IOM0) at 
Ds ~ 200h^ Mpc and the intermediate mass BH bi- 
naries of mass (IO^Mq + IO^Mq) as well as (10^1/0 -f 
IO^Mq) at Ds = 3000h"^ Mpc (Hubble distance) as the 
sources. Note here that we adopt h = 0.7. For each mass 
case we randomly distribute 10* binaries on the celestial 
sphere at Ds, and we show the probability distribution 
of the angular resolution for these sources. 
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In Fig. 2, we show the relative probabihty distribution 
of the angular resolution for the various binary masses 
Mi^2z = 1-4Mq, 10^~^Mq in the case of 1 yr observation. 
The solid lines and the dashed lines show AOs and A(f)s, 
respectively. We also show the signal to noise ratio S/N 
and the estimation error in the coalescence time Ate as- 
suming the template is accurate enough. S/N is indepen- 
dent of 6s and (f>s since the phase factor e'^'^"^^^ becomes 
unity in Eq.Q while Ate is found to depend on 9s and 
0s very weakly. In general, the errors (A6'5, A^c) 
simply scale as {S/N)~^. As shown in Fig. 2, the an- 
gular resolutions are typically ^ 0.1 — 10 arc minutes 
for the neutron star binaries and the black hole bina- 
ries. This is about 10 — 1000 times better than that of 
LISA (~ 1 deg). The Doppler modulation {(fiD oc / in 
Eq.QJ) effect is larger than the rotational modulation for 
the high frequency band like DECIGO/BBO. It is shown 
(Cutler & Vecchio 1998; Moore & Hellings 2002; Taka- 
hashi & Seto 2002), that the angular resolution scales 
as oc f~^{S/N)~^ for the monochromatic sources. Thus 
the angular resolution of DECIGO/BBO should be about 
100 times better than that of LISA if the S/N is equal. 
This explains the accuracy of ~ 0.1 — 10 arc minute. 

We comment the results for 10 yr observation. In this 
case, the angular resolution is 2 — 4 times better than the 
case for 1 yr observation in Fig. 2. 

In Table. 1, we show the S/N and the estimation er- 
rors of the chirp mass Mz, the reduced mass fXz, the 
coalescence time tc and the distance to the source Ds- 
The results are presented for the various binary masses 
(1.4, 10^~^Mq) for lyr (the upper Une) and 10 yr (the 
lower line) observation. The errors simply scale as 
{S/N)~^ for the change of the distance to the source. 
For the accuracy of the chirp mass AA4z ~ lO"'^ and 
the reduced mass A^z ~ 10"'', each mass of the binary 
is determined within 10~^. The distance to the binary 
is determined by ~ {S/N)~^. 

4. DISCUSSIONS 

The distance to gamma ray bursts (GRBs) was not 
determined before 1997. The minimum possible distance 
was 100 AU and the maximum one was the Hubble 
distance (Fishman and Meegan 1995) so that many the- 
oretical models were equally possible. The main reason 
why the distance to GRBs was not determined is that 
the accuracy of the spatial position of GRBs using the 
gamma ray observation was at most ^ deg so that many 
possible host galaxies exist in an error box of GRBs. To 

* The situation for GRBs is a little different from that for coa- 
lescing compact binaries. For compact binaries, we can determine 
the distance without identification of a host galaxy. But there are 
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find the host galaxy of GRB was akin to searching a nee- 
dle in a haystack. In 1997 Italian-Dutch satellite Beppo- 
Sax succeeded in obtaining the high resolution (~ 2 arc 
minutes) X-ray images of GRBs which led to the determi- 
nation of redshifts and host galaxies of GRBs (Costa et 
al 1997). This clearly demonstrates that the accuracy of 
an arc minute or so is indispensable to study the objects 
in full details. 

In this Letter we showed that spatial position of co- 
alescing binary neutron stars or black holes of mass 
~ IOMq at the distance of ~ 300 Mpc can be determined 
within an arc minute from one year observation. This ac- 
curacy is comparable or better than that of Beppo-SAX 
so that we may expect the rapid progress of research like 
in GRBs after 1997. ^ In Beppo-SAX case the accurate 
position of GRBs was obtained a few hours after the GRB 
event. In our case the observation period is 1 — 10 year 
before the final merging to black hole so that we would 
distribute the alert of the event beforehand. 

Let us consider the observation from 1 yr before the 
merging up to a day before the merging. In Fig. 4, we 
show the angular resolution A(f>s for the neutron star bi- 
nary (1.4 + I.4M0) at 200/i"i Mpc at each time. We 
integrate the waveform in Eq.(|3J from 1 yr before the 
final merging to three months (the dotted line), a month 
(the short dashed line), a week (the long dashed line) 
and a day (the dot-dashed line) before the merging. The 
solid line shows for the case of full one year observation 
up to the cut-off frequency. We see that the accuracy is 
~ 10 arc minutes at three months before the final merg- 
ing while ~ 0.4 arc minute at a day before so that the 
accuracy of the position depends on the time before the 
final merging. 

In conclusion DECIGO/BBO in ^ 2020 can determine 
the angular position of the neutron star binary and the 
black hole binary at ^ 300 Mpc within an arc minute 
before the final merging to the black hole. In a sense 
DECIGO/BBO will truly open gravitational wave as- 
tronomy since the positional accuracy is comparable to 
X-ray telescope on Beppo-SAX. 
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Fig. 1. — The strain sensitivity of DECIGO/BBO. We are using a smooth curve and neglecting the wavy behavior of the transfer 
function for simphcity. 
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Table 1 

s/n and the estimation errors of the chirp mass mz, the reduced mass hz, the coalescence time tc and the 

DISTANCE TO THE SOURCE Ds- THE RESULTS ARE PRESENTED FOR THE VARIOUS BINARY MASSES (1.4, 1O^-3M0) FOR lYR 

(the UPPER line) AND 10 YR (tHE LOWER LINE) OBSERVATION. 
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Fig. 2. — The relative probability distribution of the angular resolution of DECIGO/BBO in the case of 1 yr observation. For each 
mass case we randomly distribute 10* binaries on the celestial sphere at Ds (the Hubble parameter h = 0.7), and we show the probability 
distribution of the angular resolution for these sources. The solid lines and the dashed lines show A65 and A4>g, respectively. The signal 
to noise ratio S/N and the estimation error in the coalescence time Ate are also shown. S/N is independent of 65 and (f>s since the phase 
factor e"^o{f) becomes unity in Eq.Jl} while Ate is found to depend on 9g and 05 very weakly. 
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Fig. 3. The angular resolution A4>s for the neutron star binary (1.4 + 1.4Mq) at 200/1^''" Mpc (h = 0.7). The observational period 
from 1 yr before the final merging to three months (the dotted line), a month (the short dashed line) , a week (the long dashed line) and 
day (the dot-dashed line) before the merging. The solid line shows for the case of full one year observation up to the cut-off frequency. 



